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We investigate the possibility of unconventional spin density wave (SDW) in the attractive 
Hubbard model with finite spin polarization. We show that pairing and density fluctuations in- 
duce the transverse d-wave SDW near the half-filling. This novel SDW is related to the d-wave 
superfluidity induced by antiferromagnetic spin fluctuations, in the sense that they are con- 
nected with each other through Shiba's attraction-repulsion transformation. Our results predict 
the d-wave SDW in real systems, such as cold Fermi atom gases with population imbalance 
and compounds involving valence skipper elements. 
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It is well known that the Pauli's paramagnetic ef- 
fect (PPE) strongly affects superconductivity. Under a 
strong magnetic field where the conventional BCS state 
no longer exists, 1,2 the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state has been predicted to occur, 3,4 where 
the Cooper pairs have a finite center-of-mass momen- 
tum. The FFLO state has been recently reported in 
heavy fermion compound CeCoIns 5-7 and organic su- 
perconductors such as k-(BEDT-TTF) 2 Cu(NCS) 2 8 and 
A-(BETS) 2 FeCl 4 . 9 

Recently, the PPE has been also discussed in a super- 
fluid 6 Li Fermi gas, where spin polarization is realized 
by trapping atoms in two hyperfine states with differ- 
ent populations. 10 ' 11 In a cold Fermi gas, we can tune 
various physical parameters, such as pairing interaction 
and particle density. In addition, loading a Fermi gas on 
an optical lattice produced by the standing wave of laser 
light, 12,13 one can realize the attractive Hubbard model. 
Thus, cold Fermi gases are useful for the study of PPE in 
the wide parameter region of both continuum and lattice 
systems. 

The attractive Hubbard model (with a nearest- 
neighbor repulsion) is also applicable to compounds 
involving the valence skipper elements, such as 
Bai_ K Kj;Bi03. 14-16 Indeed, the phase diagram in the 
x — T plane of Bai-^K^BiOa has been nicely explained 
based on this model. Superconductivity in these mate- 
rials appears by carrier doping into the charge density 
wave (CDW) state. 17 

In this paper, we investigate the PPE in the attractive 
Hubbard model. When the spin polarization is finite, we 
predict the d-wave spin density wave (ciSDW) near the 
half-filling, where the SDW order parameter has a d-wave 
symmetry. To show this, we employ two approaches. The 
first approach uses Shiba's attraction-repulsion transfor- 
mation. 18 In the second approach, we numerically deter- 
mine the stable region of dSDW within the fluctuation 
exchange (FLEX) approximation. 19 



We consider the attractive Hubbard model under an 
external magnetic field h in the ^-direction, given by 



H 



[ c L c j>+ 1i - c - 



ha z 



(1) 



Here, c\ a is the creation operator of a fermion with spin- 
<r, m a = Cj^Cia, and a z = ±1 as a =f, |. t describes 
the nearest-neighbor hopping, and the summation in the 
first term is taken over nearest-neighbor pairs. U > 
is the on-site pairing interaction, and fj, is the chemical 
potential. In the case of two-component Fermi atom gas, 
jit = (/if + /tjJ/2 and h = (/if — /*j.)/2, where /i CT is the 
chemical potential of atoms with (pseudo)spin-cr. 

To discuss the possibility of dSDW, we assume a bipar- 
tite lattice (where the lattice points can be divided into 
two sublattices). In this case, we can always transform 
(1) into the repulsive Hubbard model by the particle- hole 
transformation 18 



ST ' "if u *i ^ • (2) 

where Q is chosen so as to satisfy exp(iQ • R) = — I for 
any translation vector R between two sublattice sites. In 
particular, Q = (tt/cl, tt/o) in the case of square lattice 
(where a is the lattice constant). The resulting repulsive 
Hubbard model is given by 



H = —t 



h.c. 



U 
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U 



(3) 



One finds from (1) and (3) that roles of [i and h are 
exchanged between the two. This means that the mag- 
netic (doping) phase diagram in the attractive Hubbard 
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Tabic I. Correspondence of order parameters between attractive and repulsive Hubbard models. CDW: charge density wave. T-SDW: 
transverse spin density wave. L-SDW: longitudinal spin density wave. SF: supcrfluid phase. 7 P = cos(p^a) — cos(p y a) is the basis 
function of d x 2 _ y 2-wave symmetry, is a coupling constant which induces d-wave SDW or d-wave SF. In this table, q is an arbitrary 
momentum. In the 2D square lattice, we take Q = (it / a, ir / a) . 



attractive Hubbard model 


repulsive Hubbard model 
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model is mapped onto the doping (magnetic) phase dia- 
gram in the repulsive Hubbard model. In this mapping, 
the order parameter in each phase is also transformed by 
the particle-hole transformation in (2). We show some 
examples in Table I. 

Since the discovery of high-T c cuprates, it has been 
widely recognized that the ((quasi) two-dimensional) re- 
pulsive Hubbard model has the (4 2 ^2 -wave supercon- 
ducting phase near the half-filling. This originates from 
a pairing interaction mediated by antiferromagnetic spin 
fluctuations. 19-24 Using this, we find from Table I that 
the transverse c?SDW phase appears in the attractive 
Hubbard model when spin polarization is finite. The pos- 
sibility of rfSDW based on the attraction-repulsion trans- 
formation have been pointed out by Ho et al. 25 

The possibility of dSDW has been discussed in the 
context of hidden ordered state in URu 2 Si2. 26 In the 
mean field theory, the ciSDW is known to be triggered 
by a nearest-neighbor repulsion. In the attractive Hub- 
bard model, this interaction is considered to be induced 
by higher order processes mediated by pairing and den- 
sity fluctuations. 

The (ZSDW state has various interesting properties. 
The (transverse) (iSDW in the 2D square lattice is char- 
acterized by the complex order parameter <f>p Q+q = 

5 , d7pSk( c kT c k-Q-ql)7k, where 7 P = cos(p x a)-cos(p y a) 
and gd is the interaction for inducing dSDW. Since 
Q+q nas n °des, the low-temperature specific heat be- 
haves as C cx T 2 . Because of the momentum dependent 
order parameter, the rfSDW state has no magnetic or- 
dering. Instead, it is accompanied by alternating circular 
spin current, 2 given by 
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where a l a/3 is the Pauli matrices, and I — x,y,z describes 
the spin component flowing from the j-th site to z-th site. 
In the 2D square lattice, the (iSDW state with q = gives 



proximation 19 is known to reasonably describe the d- 
wave superconducting phase transition temperature T c . 
In this paper, we apply this approximation to the attrac- 
tive case in (I), to examine the stability of the cZSDW 
state. The FLEX approximation has been recently used 
to determine T c in the BCS-BEC crossover regime of the 
attractive Hubbard model without spin polarization. 28 

In the FLEX approximation, the self-energy 
S,j(k, iu) n ) in the single-particle Green's function 
Gj 1 (k,k*; m ) = iu m - e k + p, + ha z - £ CT (k,iw m ) (where 
£k = —2t[cos(k x a) + cos(k y a)]) is given by 

£ CT (k,iw ro ) - [T (pp) (q,w n )Ga(q- k,w n - iuj m ) 

q.n 

+r^ ) (q, iv n )G a (q + k, \v n + \uj m ) 

+T < f h) (q, w n )G s (q + k, w n + iw m )] . (6) 

Here, uo m and v n are fcrmion and boson Matsubara fre- 
quencies, respectively, and {f , 1} = {I, T}- and 
T^ ph ^ (j = 1,2), respectively, describe the vertex func- 
tions in the particle-particle and particle-hole channels, 
having the forms 

r(pp)(q>„) = =g- -, (7) 

r ( P h), . , c/ 4 n<P h) (q,i^)ni p g h)2 (q,i^) 

IV (q,i^) - - — — (T^y— "3— -. I 8 ) 
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The diagrammatic structures of (7)- (9) are shown in 
Fig.l The polarization functions in (7)-(9) are given by 

(q, w n ) = T G a (k, iu m )G a , (q - k, w n - iluJXO) 

k.m 



(Jfj+iJfj) = it<f>daim(6 Xt , Xj ±a5 Vi , Vj - 5 Xt , Xj 5 Vi , Vj ±a) Cfq.'^) - -r^G^k^^tq + k.^ + ^l) 



x exp[iQ • Ri 



(5) 



Here, Ri = (a;,,t/j), and we have simply written the or- 
der parameter as 0p S qJo — Qdlp^dSDW- We find from (5) 
that the argument of 4>dSDW determines the spin compo- 
nent of the spin current. 

In the second approach, we determine the phase dia- 
gram of the attractive Hubbard model in the 2D square 
lattice. For the repulsive Hubbard model, the FLEX ap- 



k,m 

We briefly note that only r( pp ) involves the first and 
second order terms with respect to U to avoid double 
counting. We calculate the particle density for each spin 
n a from the condition n a — T^k m ^(k, io> m )e~ 1Wm0 . 

We determine the superfluid transition temperature 
T c when the condition 5 = 1- C/n^ p) (q,0) = 0.001 is 
achieved. Although the criterion for T c should be S = 0, 
we take the finite value of 5, because the convergence 
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Vertex functions of particle-particle channel r' pp ' and 
particle-hole channel r^ p ^, r2 Ph \ and rl ph \ which are involved 
in the self-energy (6) and SDW Eliashbcrg equation (12). The 
solid and dashed lines describe G a and U, respectively. 



of (6) is not guaranteed when S — in the 2D sys- 
tem. We note that the superfluid (or FFLO) instability 
and CDW (or incommensurate CDW) instability occur 
simultaneously at the half-filling, because of the sym- 
metry property of the model Hamiltonian in (1). Since 
our formalism treats pairing and density fluctuations in 
a consistent manner, it satisfies this required condition. 



Indeed, one finds that II 



(pp) 



nl p t h) (Q + q ,o) 



n[P h) (Q + q, 0) at the half-filling. 

The SDW phase transition temperature is de- 
termined by solving the linearized Eliashbcrg 
equation in terms of the SDW order parame- 
ter. It is obtained from the anomalous self-energy 
in the Green's function G^i (k, k + Q, w m ) = 

- Io /T dr(T T (cl (r)c k+Qi (0)))e i — , as 

, iu m ; k' , iu> m > ) G-f (k' , iuw ) 

k' ,'tn' 



xG|(k' + Q, iw m ')0(k, iu)„ 



(12) 



Here, we have included both k- and w m -dependence of 
(transverse) SDW order parameter <f>(k, iu m )- The inter- 
action V(k, ie m ; k', ie m >) in the ciSDW channel has the 
form 



;k', 



r< ph) (q>„) 



= r(pp)(k' + k + Q,i^ + i Wm ) 

+4 ph) (k'-k> m / -iw m ),(13) 

-tf 3 n<f ) (q> n )n<f ) (q,ii/ n ) 



-(14) 



" 3 ™ 1 - (q, i^)nj ph) (q, w n ) 

The SDW phase transition is achieved when the max- 
imum eigenvalue (= A max ) in (12) reaches unity. To 
seek out A max , we use a power method. 19 The result- 
ing (/>(k, iuj m ) at the half-filling has the <i-wave symmetry 
cos(k x a) — cos(k y a) 1 and is an even function with respect 
to oj m . In this regard, the odd-frequency SDW is not 
obtained unless the odd-frequency pairing is stabilized 
in the repulsive Hubbard model in (3). Thus, when one 
considers a system obtained from the extended repulsive 
Hubbard model in which the odd-frequency pairing ex- 
ists 29,30 by the attraction-repulsion transformation, an 
odd-frequency SDW is expected to be realized. 
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Fig. 2. Calculated phase diagram of the attractive Hubbard 
model in the 2D square lattice. The solid line with filled circle 
and dash-dotted line show T^SDW and T c , respectively. At the 
half-filling, s-wave superfluidity and CDW are degenerate below 
T c even when M = n-f — rij is finite. Under the condition such as 
U = At and n = 1, FFLO and Incommensurate CDW transition 
temperature is lower than T c for q = and T^SDW _ 
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Fig. 3. Calculated dSDW phase transition temperature T c SDW 
and superfluid transition temperature T c as a function of car- 
rier doping. In this parameter region, FFLO phase appears with 
center-of-mass momentum q = (±37r/16a, 0), (0, ±37r/16a) at 
the T c . 



We briefly note that, since the spin rotational sym- 
metry is broken when h ^ 0, the longitudinal <i-wave 
spin susceptibility does not diverge, in contrast to the 
transverse one. As a result, the longitudinal g?SDW is 
not obtained in the present model. 31 

The cZSDW vector Q deviates from Q = (ir/a,ir/a), 
when the carrier density is away from the half-filling. In 
addition, even-frequency SDW and odd-frequency SDW 
may mix with each other, because the non-diagnal ma- 
trix elements between even and odd frequency gap re- 
main in the interaction V(\t, iu) m ; k', iu) m i) due to the 
particle-hole asymmetry of the density of states. How- 
ever, in this paper, we only consider the (iSDW order 
with Q = (n/a, ir/a) and even frequency type for sim- 
plicity. In numerical calculations, we take 32 x 32 lattice 
points and retain the fermion and boson Matsubara fre- 
quencies tO n max = TOmax = 1024. 
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Fig. 4. Calculated dSDW phase transition temperature T c dSDW 
as a function of pairing interaction U. 



Figure 2 shows the phase diagram of the attractive 
Hubbard model in the 2D square lattice. When the spin- 
polarization M = nj — nj is small, the s-wave superfiuid 
phase appears. This phase corresponds to the antiferro- 
magnetic phase in the repulsive Hubbard model. In the 
vicinity of this superfiuid phase, we find the rfSDW phase, 
corresponding to the <i-wave superconducting phase in 
the repulsive case. The maximum dSDW phase transi- 
tion temperature T c dSDW equals 0.025i. 

Fig. 3 shows the doping dependence of T c dSDW when 
M = 0.15. The dSDW disappears at \n - 1| ~ 0.05, but 
q = superfluidity or FFLO superfluidity appears, de- 
pending on the M. This implies that the rfSDW is not so 
robust against the carrier doping as the superfiuid phase. 
Therefore, the symmetric density of states with respect 
to the Fermi level is favorable for the dSDW state. This 
dSDW-SF transition is discussed in the context of rfSC- 
AF transition induced by the magnetic field in the repul- 
sive Hubbard model. 32 We also note that the strong pair- 
ing interaction is crucial for the dSDW state, as shown 
in Fig. 4 (although j' c dSDW again decreases when the in- 
teraction is very strong (U S> 4i)). 

To conclude, we have discussed the possibility of 
dSDW state in the attractive Hubbard model with spin 
polarization. Then, we have predicted this novel SDW 
state based on the attraction-repulsion transformation, 
as well as numerical calculation within the FLEX ap- 
proximation. 

A cold Fermi gas loaded on an optical lattice is a strong 
candidate for the system to realize dSDW. In this paper, 
we considered a two-dimensional system in determining 
the phase diagram. From the previous work 33,34 on the 
quasi-two-dimensional repulsive Hubbard model, one ex- 
pects that the three-dimensionality suppresses the dSDW 
and hides it in the FFLO phase. Thus, the dimensional- 
ity of system should be close to two dimension for our 
purpose. In this regard, the optical lattice is very con- 
venient, because the dimensionality can be tunable. For 
the observation of the dSDW state, the photoemission 
spectroscopy 35 would be useful, where the single-particle 
spectral weight is expected to show the d x 2_ y 2-w&ve gap 



structure. 

The valence skipping compounds, Bai_ x K x Bi03 17 and 
BaPbi„ x Bi x 03, are also other candidates. Their phase 
diagrams are well described by the extended attractive 
Hubbard model involving the nearest neighbor repul- 
sion. In the mean-field theory, this additional repulsion is 
known to be favorable for the stabilization of the c?SDW 
state. Thus, it is a possible scenario that the rfSDW ap- 
pears in the vicinity of the CDW or superconducting 
phase under a magnetic field. 

We note that the stability of dSDW against the or- 
bital effect (which always exists in charged systems) is 
still unclear, although the d-wave density wave state is 
known to be insensitive to this. 36 We also note that the 
dSDW may coexist with superfiuid state or CDW near 
the phase boundary below T^ SDW . These are important 
future problems to understand the physics of unconven- 
tional spin density wave. 
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